Exploration and Trapping of Mortal Random Walkers 



cn 

o 

(N 

< 

oo 
(N 






S. B. Yustc^, E. Abad^, and Katja Lindcnbcrg'^ 

^^' Departamento de Fisica, Universidad de Extremadura, E-06071 Badajoz, Spain 

^^' Departamento de Fisica Aplicada, Centra Universitario de Merida, 

Universidad de Extremadura, E-06800 Merida, Spain 

^^' Department of Chemistry and Biochemistry, and BioCircuits Institute, 

University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0340, USA 

(Dated: April 30, 2013) 

Exploration and trapping properties of random walkers that may evanesce at any time as they 
walk have seen very little treatment in the literature, and yet a finite lifetime is a frequent occurrence, 
and its effects on a number of random walk properties may be profound. For instance, whereas the 
average number of distinct sites visited by an immortal walker grows with time without bound, 
that of a mortal walker may, depending on dimensionality and rate of evanescence, remain finite or 
keep growing with the passage of time. This number can in turn be used to calculate other classic 
quantities such as the survival probability of a target surrounded by diffusing traps. If the traps 
are immortal, the survival probability will vanish with increasing time. However, if the traps are 
evanescent, the target may be spared a certain death. We analytically calculate a number of basic 
and broadly used quantities for evanescent random walkers. 

PACS numbers: 05.40.-a,05.40.Fb,02.50.-r,46.65.+g 
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Random walk models provide a quintessential ap- 
proach to transport and related processes in condensed 
media. Random walks have been studied for more than 
a century - indeed, as an antecedent, there is reference to 
probability and statistical inference in biblical texts [l|. 
It is therefore surprising to find important problems in 
this arena that have not yet been explored, especially 
ones that are broadly applicable and that can be dealt 
with analytically. This letter deals with one class of such 
problems associated with the territory explored by mortal 
or evanescent random walkers (we use the two terms as 
synonyms). Mortal walkers or mortal diffusing particles 
may disappear in the course of their motion. This disap- 
pearance may, for instance, be the result of a finite walker 
Hfetime such as in a unimolecular reaction or a natural 
decay process. Other examples of disappearance events 
may arise from an encounter of a walker with another 
walker leading to the annihilation of one or both, as may 
occur in a chemical reaction such as radical recombina- 
tion, in exciton trapping in photosynthesis, or in growth 
by aggregation. Many references to suchphenomena can 
be found in [24j|, and most recently in [5|. 

The statistical properties of the territory explored by 
immortal random walkers as a function of time have been 
studied in a variety of contexts |6l4ld|. So have related 
quantities such as the probability of return to a given lo- 
cation. In turn, there are further connections between 
these and reaction kinetic quantities such as the sur- 
vival probability of a target particle surrounded by diffus- 
ing traps ( "target problem" ) [ll| - il4J | . Early last century 
these problems experienced an enormous surge in the lit- 
erature with the pioneering work of Smoluchowski in the 
context of diffusion-limited chemical kinetics. In more 
recent times, a resurgence of interest in random walks 
started with the classic works of Scher et al. on stochastic 



transport in amorphous condensed media |15l . Il6| | . Traps 
or defects in these media are slowed down by the disor- 
dered environment and consequently experience so-called 
anomalous diffusion. This slowing down, when incorpo- 
rated in the "defect diffusion model," leads to stretched 
exponential relaxatioii, which turns out to be widely 
ubiquitous in nature [51, [iTl |l8[. The periodic surge of 
extraordinary interest has once again been proved by the 
plethora of recent books and chapters on anomalous dif- 
fusion models as a descriptive tool of crowded disordered 
condensed systems 19|, l20[ • 

Statement of the problem. We consider a sym- 
metric nearest-neighbor random walk on a rf-dimensional 
lattice, that is, a Polya walk. The walker steps at dis- 
crete times t„. where n is the number of steps. We will 
also consider the continuous version of the problem, a 
diffusive process in a continuous medium taking place in 
continuous time. Our first goal is to calculate a quantity 
which can then be used to calculate many others: the 
average number of distinct sites visited by an evanescent 
walker up to time i„, S"* (the walker may die at an earlier 
time, in which case S"* stops growing). The correspond- 
ing continuum quantity is the average volume Vf of the 
Wiener sausage generated by mortal particles up to time 
t. The asterisks denote evanescent particles (the corre- 
sponding quantities for immortal walkers are indicated 
without an asterisk). In turn, these results can be used 
to address other classic problems, now for mortal walk- 
ers. Perhaps one of the most interesting arises from the 
well-known connection between the survival probability 
up to a given step number or time of a target particle sur- 
rounded by a concentration of diffusive evanescent traps. 
For immortal walkers, in the discrete problem this sur- 
vival probability is (/)„ = exp[— p(S'„ — 1)] and in the 
continuous case it is (f>{t) = cxp{—cR'^Vt), where p and c 



denote the density of walkers in appropriate units and R 
is the radius of the target (assuming point traps; other- 
wise R is the sum of the radii of the target and a trap). 
These relations persist for mortal walkers. 

We start with the number of distinct sites visited by 
an evanescent walker. To arrive at this quantity we intro- 
duce the probability P^ n-mi^\^') o^ finding an evanes- 
cent walker at lattice site s after having taken n — m 
steps if the walk started at site s' at step m. The initial 
condition is Pg*„(s|so) = Ssso- This probability should 
be thought of as the outcome of carrying out the ex- 
periment repeatedly, starting the walker at site sq and 
counting the fraction of realizations that arrive at site s 
at step n. Alternately and equivalently, if a number of 
noninteracting walkers all start at step at site sq, this 
is the fraction that arrive at site s after n steps. This 
probability is related to the corresponding well-studied 
probability Pm.nis\s') for immortal walkers via the rela- 
tion P^,„(s|s') = [pin)/p{m)]Pm,n{s\s'), where p{n) is 
the fraction of realizations for which the walker has not 
evanesced up to step n or, alternately, the concentration 
of walkers that have not evanesced up to that step (with 
p{0) = 1). We also introduce F*{s\so), the probability 
that the evanescent walker arrives at site s for the first 
time at step n if the walker started at site sq at step 
n = 0. The probabilities P^ „ and F* are related in the 
same way as for immortal walkers: 

n 

F;o(s|so) = Sss„Sno + Y.F*{s\so)Pl,,_^{s\s), n > 0. 
i=i 

(1) 
Let A* denote the average number of new sites (sites 
never visited before) visited by the n*'' step of an evanes- 
cent walk, with AJ = 1. Then 



5: = 5]A* = 5]p(j-)A,. 

j=o J=0 



(2) 



We define the generating function of any n-dependent 
quantity A„(-) by A(-; ^ = EZo ^"(OC"- The generat- 
ing functions of S** and A* arc then related by 



s*{0 



A*(0 



(3) 



On the other hand. A; = J^s^so K{s\so), n^O. Mul- 
tiplying by ^", summing over n, and reversing the order 
of summation yields 



A*{0 = 1+J2f*{s\s„;0. 

s^so 



(4) 



In order to go further we need to specify particu- 
lar forms of evanescence. We consider exponential and 
power-law decay of the concentration of walkers. The 
former is the typical unimolecular decay that describes 



spontaneous death; the latter is typical of more complex 
chemical reactions [J, [21 1 . 

Exponential evanescence. With exponential 

evanescence, p{n) = exp(— An). In a walk of immortal 
walkers on a regular lattice Pj,„(s|so) = Pn-j(s|so), that 
is, the walk is time invariant. Exponential evanescence is 
the only form of evanescence that preserves this property 
for P* . It then follows from Eq. (P that 



F*is\so;0 



P*is\so;0-6s 
P*(0;e) 



(5) 



where the translational invariance of our lattices implies 
that P*(s|s;0 = P*(so|so;C) = P*(0;0- Then, from 
Eq. (O with Eqs. © and (g]) one finds 



,.(,) 1 ^Pn.|.o;0 



P*(0;e) 



(6) 



and then, using the abbreviated notation ^ 



,-A 



L 



(1 - 0S*{0 = [(1 - O^(0;O]"^- Here we have used 
the relations J2 s P* is\so; = p{0 = 1/(1 - I) and 
P*(0;C) = P(0;0- Lattice Green functions P{0;O are 
well known for the most relevant d-dimensional lattices 
[111 |22| , which then allows us to find a number of results 
for the evanescent walk. One of the interesting quantities 
obtained from the expansion of S*{^) in a power series 
is the average number of sites visited up to time t„ by a 
mortal Polya walker with exponential evanescence. For 
immortal walkers 5^ ^^ oo as n — s- cxd. For mortal walk- 
ers it is finite in the case of exponential evanescence, and 
is given by 



1 



1 



1-e 



-A 



P(0;e-^) 



(7) 



Specific values of S^ depend on dimension and type of 
lattice. In dimension d = 1, P(0;O = (1 - C^)"^^^ 



so that S* 



[(1 



^)/(l 



')] 



1/2 



For a two- 



dimensional square lattice P(0;^) = 2K{^)/tt, where 
A'(-) is the elliptic integral of the first kind. 

From the known asymptotic behaviors of P(0;^) as 
^ ^- 1~ and the fact that ^ = e^'^^, one can arrive at 
the large-n behavior of S** for A — !■ 0, that is, for slow 
evanescence. Focusing first on the leading asymptotic 
contribution, we note that for d = 2, P(0;^ -^ 1~) ~ 
A/7rln[P/(l — ^)], where the constants A and B depend 
on the type of lattice [uj . This behavior in Eq. ([7]) then 
yields S^ ~ 7r/[A Alog(B/A)] as A ^- 0. In dimension 
d > 3 the probability that a walker returns to the origin 
is given by 7^ = 1 - 1/P(0; 1). Hence, 5^ - (1 - 7^)A-l 
as A ^^ 0. 

The way in which 5* approaches S^ for large n can 
be obtained from the subdominant behavior of P(0;^) 
as ^ ^' 1~. For three-dimensional lattices P(0;^) = 

E™=o(-ir"™(i - 0™/' = [E™=o^™(i-0"/']"\ 
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FIG. 1. S^ vs n for the simple cubic lattice and, from top 
to bottom, A = 0, 0.01, 0.02, 0.05, 0.1, 0.2. Solid lines: Eq. P; 
broken lines: S^ given by Eq. ((7|; squares: exact values ob- 
tained by identifying the 100 first coefficients in the ^-power 
expansion of S*{^); circles: simulation values for 10''' runs. 
The excellent performance of the asymptotic expression ^ 
even for small n is shown in the inset. For the simple cubic 
lattice Mo = 1/(1 - 7^) ~ 1.51639 and m = 3'^^'^ /{ttV2). 
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FIG. 2. Sxi'~ vs S'^ and A for several values of r and A and 
three different lattices. Symbols: numerical simulations for 
d = 1 (triangles), d = 2 (square lattice, circles), and d = 3 
(cubic lattice, squares) for 10^ runs. The values of S^^ from 
the simulation of S„ with n sufficiently large to observe no 
change in at least three significant figures. From left to right: 
A = 0.1,0.05,0.03,0.01,0.05,0.001, with from top to bottom 
r = 2, 3, 4, 5. The straight lines of slope (r — 1) through the 
origin are the theoretical predictions. 



where the Urr, and Vm are known for a number of lattices 



111 |23| . Using the second expression in the result follow- 
ing Eq. ^, expanding in powers of ^, and retaining only 
the first two terms leads to 



o„ 



1 e-^("+i) Wi4-a(?i+ 1,1/2) 



Uq 1 



= -A 



(8) 



for n — > oo , where /^(a, b) is the regularized Beta func- 
tion. This asymptotic expression turns out to be surpris- 
ingly accurate even for relatively small n and for A's that 
need not be extremely small (see Fig. [1]). In fact, the 
results for A close to zero are so good that it is possi- 
ble to find the large-n asymptotic expression for Sn by 
taking the limit A — ?► of Eq. ([8]). Expanding this re- 
sult in powers of n yields a series whose first three terms 
(proportional to n, n^/^, and n^) are identical to those 
obtained by expanding the exact result for Sn [23|. Dif- 
ferences only appear in the fourth term, proportional to 
n-i/2. 

A quantity related to S** is Sn , the average number 
of sites revisited at least r times by an evanescent walker 
in an n-step walk. Following the procedure of MontroU 
and Weiss [23[ , it is not difficult to find that the generat- 
ing function for this number for exponentially evanescent 
mortal walkers is given by 



S*ir)(^^-J 



1 



PiO;0 



r-l 



S*iO- 



(9) 



From here one finds formulas for Sn 
[23|. For instance, in dimension d 



(r) 



in terms of S'* 
1 one obtains 



5*^'^ = ^; - 1 - e-^, Sn^^^ = 25; - e-2A5* _ i _ 2e-^, 
etc. The average number of sites visited r times before 
the walker dies is in any dimension given by So6 = 
[(1 - e-^)^^;,]'^-^ {S^y. This result is exhibited as a 
function of r and of A in Fig. [21 where it is also com- 
pared extremely favorably with simulation results. For 
the average number of revisits to the origin after n steps, 
fin: one finds the generating function fj,*{^) = P(0; £,) — 1. 
(We follow the convention of not counting the initial oc- 
cupancy of the origin as the first revisitation [ll|). We 
find that the average number of revisits to the origin in 
any dimension is given by /x^ = [(1 — e~^)S^]~'^ — 1, 
and the average number of visits to a site s other than 
the origin is given by /x^(s|so) = -P(s|so; e""^), one of the 
few previously known results for exponentially evanescent 
walkers (see Sec. 3.2.4 of Ref. [llj). 

Power law evanescence. Power law evanescence of 
random walkers is given by p{n) = (l + Art)^^ with A > 
and /? > 0. Here it is convenient to directly use the 
relation ([2]) and rely on the knowledge of A^, for large 
and small n for the most common lattices [llj, |23[ . For 
example, A„ ~ (1 — 7^) A 



^;3„-/3 



(l-hCn-i/2 + ---) for 
three-dimensional lattices and large n. Because p(n) ~ 
{\n)~^ for large n, one sees immediately that S'^ is finite 
for (3 > 1. For slow evanescence (A -^ 0) we find 






1-n 

(^^T)A' 



/3> 1. 



(10) 



For /3 < 1 the result for slow evanescence is quite differ- 



ent. For large n we find 



X-^n^-^, 0</3<l. 



(11) 



For the marginal case /3 = 1 one gets 5* ~ (1 — 
7^)A~^ log n. For (3 = (no evanescence) one recovers the 
classical result whereby Sn is proportional to the number 
of steps n |ll| . We thus see that the average number of 
distinct sites visited by a mortal walker before it dies is fi- 
nite for /3 > 1 , whereas this quantity is infinite for /3 < 1 . 
This is true for d-dimcnsional lattices with d > 2. How- 
ever, for the one-dimensional lattice the critical value is 
/3-1/2. 

Mortal Brownian particles in continuous space, 
and stretched exponential relaxation. It is well 
known that the results for Sn for immortal walkers can 
be used to calculate the average volume Vt of the Wiener 
sausage generated up to time t by an immortal Brownian 
particle in a continuous medium. Since the relation be- 
tween Sn and Vt is purely geometric, it can immediately 
be translated to the case of mortal walkers. Therefore, 
our results for 5* can be used to find the average volume 
v^ of the Wiener sausage generated by a mortal Brown- 
ian particle up to time t. Explicitly, 5** with n ^ 1 for a 
walker in a d-dimensional simple cubic lattice with lattice 
constant £, and the Wiener sausage volume v^ generated 
by a spherical diffusing particle of radius _R ;3> ^ up to 
time t = ne'^/{2dD) are related by v^ - ^5'*/i? for d = 1 
and by v^ ^ 7d(^/i?)^S'* for d > 2, 7^ being a constant 
that depends on dimension |24| . 

This connection in turn greatly expands the interest- 
ing world of stretched exponential relaxation discussed in 
the literature for several decades. It is well-known that 
the evaluation of the number of distinct sites visited (or 
the volume explored) up to a given time is tantamount to 
the evaluation of the survival probability (j){t) up to that 
time of a fixed target particle of radius R surrounded by 
a concentration of diffusing point traps (target problem). 
The connection between the two is (j){t) ~ exp(— ci?'^wt). 
This expression also holds for evanescent traps with the 
replacement of Vt by v^ . The identification of these traps 
as defects (i.e., carriers of free volume) is the basis of 
the defect diffusion model to explain the stretched ex- 
ponential (or Kohlrausch- Williams- Watts, famously ab- 
breviated as KWW) relaxation, in which In (/>(<) ^ t^. 
However, only the values 9 = 1/2 and 6 = 1 are possi- 
ble for normal non-evanescent diffusive defects because 
Sn oc V( oc t^/^ ioT d = 1 and (and then 6 = 1/2) and 
Sn (x Vt (X t for d > 2 (and then = 1). This limited 
model [25[ was extended by Shlesinger and MontroU [26 1 
by assuming that the movement of the defects can be 
described by a CTRW model with a power-law waiting 
time V^i) ~ i"^"'^, < 7 < 1, which leads to 6* = 7/2 
for d = 1 and = 7 for d 7^ 1 [17|. That is, stretched 
exponential relaxation with 6 ^ 1/2 in this scenario is 
explained by assuming that the diffusion of the defects is 



anomalous, with anomalous diffusion exponent 7 (lead- 
ing to subdiffusion when < 7 < 1). 

The results obtained in this Letter provide another 
route for explaining stretched exponential relaxation 
even for the case of normal defect diffusion by allowing 
the defects to disappear during the relaxation process 
2ll |27| . As we have shown, different kinds of evanes- 
cence lead to different laws of relaxation. For example, 
from Eq. ^ we see that for d > 3, < oc t^^'^ for ^ < 1, 
so that one can arrive at stretched exponential relax- 
ation with exponent = 1 — (3 when the concentration of 
defects decays as a power law. Moreover, if the concen- 
tration of defects decays as p{t) ^ 1/t for large t, which 
corresponds to /? = 1, one finds that v^ oc Ini, which in 
turns leads to algebraic relaxation [28|-|30[. 

Trapping problem. The survival probability of the 
target in the target problem is frequently and appropri- 
ately used as a first approximation (the "Rosenstock ap- 
proximation" ) to the survival probability of the target in 
the so-called trapping problem in which the target dif- 
fuses and the traps are frozen [Uj. The results of this 
Letter can also be applied to this problem, now for the 
case of traps whose concentration decreases with time 

Conclusions. As noted in the introduction, it is 
surprising to find important solvable problems involv- 
ing simple random walks, but we appear to have done 
so in the case of random walkers that evanesce in the 
course of their motion. A number of classic problems 
such as the distinct number of sites visited as a function 
of time, or the survival probability of a target pursued by 
randomly walking traps, or any number of other quanti- 
ties, change dramatically when the walkers can die in the 
course of their motion. To mention but one or two such 
changes, we showed that the average number of distinct 
sites visited by an evanescing walker in n steps as n ^ cxo 
may be finite (depending on the speed of evanescence), 
whereas it is clearly infinite if the walkers live forever. 
Another, closely related to this, is the survival probabil- 
ity of a target in the presence of mobile traps. If the traps 
live forever then the target will eventually disappear with 
certainty; if the traps evanesce, then the target may be 
spared. 

We have also enriched the world of stretched expo- 
nential relaxation of a target as calculated using the 
defect diffusion model. When the defects live forever, 
stretched exponential behavior in the CTRW model is 
only obtained if their diffusion is anomalous. Here we 
have shown that the same stretched exponential behav- 
ior is obtained with normally diffusing defects provided 
they evanesce but do so sufficiently slowly. Indeed, al- 
though not addressed here, it may happen that the traps 
are subdiffusive and also have a finite lifetime. Else- 
where we have addressed consequences of finite lifetimes 
together with subdiffusion on the basis of fractional dif- 
fusion equations 3^ |33[. There are many other prob- 



lems for evanescent walkers related to explored territory 
(or distinct sites visited) that we have not had space to 
present here. They include the distribution of the dis- 
tinct number of sites visited (rather than just the mean) , 
walks with different waiting time distributions (CTRWs) , 
and Levy flights and walks. We are currently pursuing 
these and other related problems. 
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